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[Continued from page 142.] 

Thus it is proved that when any two entire polynomials are multiplied 
together, the cube of the cubic radius about the centre of forces in the product 
is equal to the sum of the cubes of the like radii in the two factors. Hence, 
if any number of such polynomials are multiplied together, the cube of the 
radius in the final product is equal to the sum of the cubes of the radii in all 
the factors. The cube of the radius in the h power of a polynomial is k 
times the cube of the radius for the polynomial itself. These propositions 
evidently hold true also for polynomials like (5), in which the place of z° 
is not at the first or left hand term. The coefficients in a product are not 
altered when all the exponents in either factor are increased or diminished 
by a constant quantity. 

Applying the above to the expression for a in (14), we write 

showing that the part within the parentheses may be regarded as the square 
of the quadratic radius divided by the cube of the cubic radius, either in the 
first power of the polynomial or in its expansion to the h power, as may be 
most convenient. The values of a and b may thus be expressed either by 
means of the coefficients X in the given polynomial, or by means of the ordi- 
nates y to the limiting curve. When the ?'s and y's are all positive, and 
represent probabilities, kb 2 (dxf is the square of the quadratic mean error e, 
and hb z (dx) 3 is the cube of what we will call the cubic mean inequality, 
which we denote by £. The constants in (14) will then be 

a = 2e 2 -f<3 j b = e 2 . (39) 
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The square of the q. m. error is of course positive, but the cube of the c. 
m. inequality is either -f or — , according as the + or — errors prepon- 
derate in forming it. If the facility of error is the same on both sides of 
the arith. mean or centre of gravity, so that + and — errors of equal am't 
are equally probable, £ will be zero, and a becomes infinite. If £ is nega- 
tive, a is also negative, and the position of the limiting curve is reversed, so 
that it lies on the left of the origin, or rather, the origin is now at the other 
side of the expanded polynomial. Since y in (25) is a function of ax, it will 
have the same value when a and x are both — , as when they are both -f- . 

The notation we adopted in (14) is such that in the expansion of the po- 
lynomial (5) to the infinite h power we have 

(quadratic rad.) 2 — b, (cubic rad.) 3 = 26-^-a. (40) 

We will now show that these values of the radii are deducible from the 
equation (25) of the gamma curve, and thus verify the proposition that the 
curve is a true limiting form of the expansion of the polynomial, since it pos- 
sesses properties which are known to characterize that expansion. The dis- 
tance from the origin to the centre of parallel forces, or centre of gravity of 
the masses y, will be 

The divisor here is unity, so that the distance sought is 

if** = apr^K^H = ^ww = * (4i) 

This agrees with (12), for by (14) we have 

ab = 2kb%dx-r-b s . 
The squared quadratic radius of the masses y about the centre of grav. is 

ijy-abfydx+(yy x ), 

that is, the divisor being unity as before, 

±f\x-abyydx= -p-L^ { ^ / Vr +1 e-°*d(ax)~2bfjaxr b e-<"d(ax) 

-\-aWf a (ax) M - l e-" x d(ax) 1 

= Wb) \ ^^^ ~ 26r (« 26 + 1 )+ aW ( ffl26 ) } 

= 6(a 2 6+l)— 2a 2 o 2 +a 2 6 2 - b, (42) 

a result which agrees with (40). Likewise for the cube of the cubic radius, 
omitting the divisor unity, we have 
1 C(x-ab\*vdx- 1 / na*b+Z)__3blWb+2) + 3ab*r\a*b+l) 1 



—163— 

= ft-\(a?b+l)(a?b+2) — Zab%a?b+l)+3a 3 b 3 — a?b 3 = ^-. (43) 

This also agrees with (40), so that the curve (25) does exactly represent 
the form of the series of coefficients in the expansion of the polynomial (5) 
to an infinitely high power, so far as the quadratic and cubic radii about the 
centre of forces are concerned. 

As shown in my former articles, a binomial p -f- q or p -f- qz, in which p 
+ q = 1 and the coefficients p and q are separated by the interval dx, has 
its centre of gravity at the distance qdx from the first term p, and the sq'd 
quadratic radius about that centre is 

(qp 2 +pq 2 ) {dx) 2 = pq(dx) 2 . (44) 

If the binomial is raised to the mth power, the centre of gravity in the 
expansion will be at the distance qmdx from the first term, and the squared 
quadratic radius about that centre is 

e 2 _ pqm(dx) 2 . (45) 

The cube of the cubic radius in the first power is 

(qp 3 — pq 3 ) (dx) 3 = pq(p-q) {dx) 3 , (46) 

and in the mth power, as we have here shown, it is m times as great, or 

C 3 = pqm(p—q)(dx) s . (47) 

Hence, when m becomes infinite, the constants in the limiting curve, ac- 
cording to (39), will be 

2 
a = -. r-r, b = pqmidxY. (48) 

We found in (22) that y is a maximum when x = ab — 1 -~a, so that the 
vertex of the gamma curve is at the distance — 1-f-a from the centre of 
gravity, and by (48) 

-\ = -l(p-q)dx. (49) 

The agreement between this result and that which I found by different 
means in Analyst, Vol. VII, p. 3, shows that the vertex of the curve 
(25) accurately represents the position of the vertex in the expanded binom- 
ial, with reference to the ordinate through the centre of gravity. 

It will often be convenient to have the origin of coordinates transferred 
to the centre of gravity. Putting x-\-ab in place of a; in (25), we have 

dx /a 2 &\° 2 Yi , aV 25 - 1 _ax /«« 

y ^m^b)W) l 1+ ^) e ■ < 50) 

A known formula for r(n) is 
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and by means of this (50) is reduced to 

K = 1+ B36 + 288(o>6f - * C -' | (62) 

» - Kvkth^T' '-■ J 

This is the equation of the gamma curve referred to the centre of gravity 
of the masses y as an origin. To illustrate the uses of the curve, we will now 
employ it in computing the principal terms in the expansion of a binomial 
to a high power. 

Mortality tables show that among persons aged 40, about one per cent 
may be expected to die within a year, so that the probability of dying with- 
in a year is .01, and that of surviving a year is .99. Suppose we wish to 
find the probabilities that out of 1000 persons aged 40, the number of deaths 
within a year will be 0, 1, 2, 3 &c. These probabilities are the 1st, 2nd 
&c. terms in the expansion of the binomial 

(p+?)» = (.99 f.01) 1000 . (53) 

The common interval between consecutive terms in the expansion being 
dx, we will take this as the unit of abscissas. By (48) we have 

a = i9rar = iw> b = -"x- 01 * 1000 = h, 

and by (52) 

logy = 1,1022148+40.23282 log (l + .04949495z)---.8863153a\ (54) 

Since the distance qrndx = 10 of the centre of gravity of the expanded 
series from its first term is in this instance a whole number, it follows that 
one term of the expansion stands exactly at that centre, where x = 0, and 
the whole series of terms will be found by putting x successively equal to 
-2, -1, 0, 1, 2, 3 ^ 

The resulting values of y are given in column (2) of Table I. To show 
the degree of accuracy attained, the true values of the terms in the expan- 
sion have been computed and set in column (1), and the differences (2) — (1) 
are also given, in units of the fifth decimal place. The computed curve in- 
tersects the true one at four points. The agreement between them is pretty 
close, and would be closer if the exp'nt m were a greater number than it is. 

Eeverting now to the more general significance of the gamma curve, as 
representing the expansion of a polynomial, we will inquire what simpler 
form of curve it approximates to. Let (52) be written 
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Table I. 
x (1) (2) (2)-(l) (3) (3)-(l) x (1) (2) (2)-(l) (3) (3)-(l) 





- 1 

- 2 

- 3 

- 4 

- 5 

- 6 

- 7 

- 8 

- 9 
-10 
-11 
-12 
-13 
-14 



12574 
.12562 
.11283 
.08999 
.06274 
.03746 
01861 
.00739 
00220 
.00044 
.00004 



.12654 
.12635 
.11311 
.08966 
.06203 
.03680 
.01834 
.00748 
.00242 
.00060 
.00011 
.00002 



80 

73 

28 

-33 

-71 

-66 

-27 

9 

22 

16 

7 



.12679 
.12055 
.10360 
.08048 
.05651 
.03587 
.02058 
.01067 
.00500 
.00212 
.00081 
.00028 
.00009 
.00002 
.00001 



105 

-507 

-923 

-951 

-623 

-159 

197 

328 

280 

168 

77 



11431 
09516 
07305 
05202 
03454 
.02148 
.01256 
.00693 
00362 
,00179 
,00084 
,00038 
,00016 
,00007 
,00003 
.00001 



.11482 
.09525 
.07282 
.05167 
03422 
.02127 
.01247 
.00692 
.00365 
.00184 
.00089 
.00041 
.00018 
.00008 
.00003 
.00001 
.00001 



51 
9 
-23 
-35 
-32 
-21 

- 9 

- 1 
3 
5 
5 
3 
2 
1 


1 



.12055 
10360 
.08048 
.05651 
.03587 
.02J58 
.01067 
.00500 
.00212 
.00081 
.00028 
.00009 
.00002 
.00001 



624 
844 
743 
449 
133 

— 90 
—189 
-193 
—150 

— 98 

— 56 

— 29 

— 14 

— 6 

— 3 

— 1 



'.log'(|) = (o 2 6— l)log'(l+^) — ax 

c -^{s-j&w-m- 



-ax 



x* , (x* l\x (x* l\/a;\2 /a 2 1\/ x\ z ,„, 

= -26 + (36-Tfe-<46-2Jy ^ST"^!^) "'^ 

Since 6 is a finite area, x 2 -i-b is in general a finite number. By (14) we 

have 

te_ lbj\xdx (g6) 

ab \bj2b' y ' 

which is in general an infinitesimal. Neglecting all terms in which x-^-ah 
is a factor, (55) is reduced to 



MJ) - -w 



y = oe-**+* b . 



Restoring the value of c from (52), noticing that when k is really infinite 

we have 

tfb = 4**f-5-&i = °o, ••• K = 1, 
we get finally 

dx 



» - V&F)^™ (67) 

the equation of a common probability curve like (1). This is the same re- 
sult we would have obtained in the first place, if we had negl'd d % y in (10). 
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If instead of retaining only dy and d 2 y, we should also retain d s y, the re- 
sulting equation, if we could integrate it, would doubtless give a limiting 
curve of still more general form, of which the gamma curve is but a partic- 
ular case. Under this view, the probability curve (57) and the gamma 
curve (52) are only first and second approximations to the actual form of an 
expansion to a high power. 

It must be observed that since dx represents the common interval between 
consecutive coefficients y in the expanded series, the abscissas x correspond- 
ing to ordinates at and near the origin will have the values 

.... — dx, 0, dx, 2dx, . . . . , 
so that x 2 -±b and x^-ab in (55) will there be of the same order of magni- 
tude, and the latter cannot be neglected in comparison with the former. 
The curve is thus rendered unsymmetrical in the immediate vicinity of the 
origin, and the maximum or vertex is thrown a very little to one side of the 
centre of gravity. But if the given polynomial (5) is such as to make & 3 =0, 
then (56) makes x-i-ab absolutely null, and the terms in which it is a factor 
disappear altogether from (55), leaving (57) as the exact result, and show- 
ing that the probability curve is a special case of the gamma curve, occurr'g 
when 1-r-a = 0. This case arises at the moment when the gamma curve 
passes from the direct to the reversed position, as noticed in connection with 
(39). The X axis then becomes an asymptote to the curve, not on one side 
only, but on both sides. Again, if 6 3 , though not absolutely zero, is quite 
small in comparison with 6 2 , making £ 3 quite small in comparison with s 2 , 
a thing which often occurs, then a is a very large number, and terms con- 
taining x-r-ab may be dropped from (55) as before. In other words, if the 
asymmetry, as measured by the c. m. inequality £, is small, the gamma 
curve does not differ materially from the common probability curve, and 
the latter may be preferred to it for practical use, as being more simple. 

But when £ is so large that a is a small number, the probability curve 
will not represent the true form of the expansion of a polynomial to a finite 
power with sufficient accuracy. Take for instance that of the binomial (53), 
for which, with dx = 1 as before, and b = 9.9, we get by (57) 

logy = 1.1030924— .02193407a; 2 . (58) 

The terms of the series computed by this are entered in column (3) of 
Table I. Their sum is unity as it should be. The differences (3) — (1) 
between the computed values and the true ones are also shown. On an av- 
erage they are more than 12 times as great as the differences (2)-(l) afforded 
by the gamma curve. 

That the expansion of a polynomial to a high power, taken as a whole, 
tends to become more and more symmetrical in form the higher the power 
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is, may be inferred from the properties of the quadratic and cubic radii, 
without regard to any precise analytical expression for the limiting curve. 
While the whole length of the expanded series increases in proportion to 
the exponent k of the power, the quadratic radius of the coefficients, about 
their centre of forces, increases only as the square root of k, while the cubic 
radius increases still more slowly, being proportional to the cube root of k, 
as seen in connection with (37). 

Suppose that the coefficients A in the polynomial (5) represent the proba- 
bilities of the occurrence of the various possible true errors of an observed 
quantity, these errors being multiples of the unit of measure Ax, which may 
be taken as small as we please, while m is a whole number so large that the 
greatest error will not exceed ± mAx. In any term 1$ the coefficient ^ is 
the probability that, in a single observation, the error which occurs will be 
x = iAx. The centre of gravity of the coeffic's, regarded as the masses of 
material p'ts ranged along the imponderable axis of X, may or may not co- 
incide with the place of X Q , but wherever it is, its abscissa x lf or lever arm 
about the place of A , is the arithmetical mean of all the possible true errors 
of a single observation, each error being taken with a weight proportional 
to the probability of its occurrence. Let e and £ denote the quadratic and 
cubic radii of the masses k, about the centre of gravity. These are the same 
as the q. m. error and c. m. inequality of a single observation, if what we 
call errors are not necessarily true errors, but only deviations from the cen- 
tre of gravity or ultimate arith. mean. (Analyst, VIII, p. 141.) 

If k such observations are taken, the possible true errors of their sum will 
be the exponents, and their probabilities will be the coefficients I, in the po- 
lynomial (6), which is the expansion of (5) to the k power. The centre of 
gravity of all the coefficients I will be approximately the place of the max- 
imum coefficient, and its abscissa, or lever arm about the place of l , will 
be kx x . This lever arm is the arithmetic mean of all the possible true errors 
in the sum of k observations, the errors being weighted for probability of 
occurrence. The quadratic and cubic radii for the masses I, about the cen- 
tre of gravity, will be 

E = e x /k, Z = Qfk, (59) 

and these are the q. m. error and c. m. inequality of the sum of k observa- 
tions. 

The probability that the true error of the sum of k observations will be 
x, is the same as the probability that the true error of their arith. mean will 
be x-s-k. If we suppose the coefficients or masses I to be set closer together, 
so that the common interval between them is reduced from Ax to Ax-s-k, 
their distribution along the X axis will represent the law of facility of error 
in the mean of k observations. The limits of possible error, which were ± 
kmAx for the sum, will be reduced to ± mAx for the mean, being: the same 
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as for a single observation. The centre of g. of all the masses I will now, aa 
before, be the approximate place of the maximum, and its abscissa, or lever 
arm about the place of l , will be reduced from kx x to a? 1 , showing that if 
there is any true error in the ultimate mean, or arith. mean of all the possi- 
ble values, weighted for probabi'ty of occurrence, it is the same for the arith. 
mean of k observations, as it is for a single observation. 

Since the distance of each of the masses I from the centre of gravity is k 
times less for the arith. mean than it was for the sum, 1-t-Jc becomes a com- 
mon coefficient of all the distances which go to make up the quadratic and 
cubic radii, which are consequently k times less for the mean than they were 
for the sum. Hence by (59), the q. m. error and c. m. inequality for the 
arith. mean of k observations will be 

e = e H-i/£, Co = C-*-**- (60) 

When k is increased, the q. m. error of the mean diminishes, being inv'sely 
as \/k, while the c. m. inequality diminishes more rapidly, being inversely 
as f/k 2 . If we take 64 times as many observations, the q. m. error of the 
mean result will be one eighth as large as before, but the c. m. inequality of 
the possible errors of the mean will be only one sixteenth as large as before. 
This goes to show that as k increases, the curve of facility of error in the 
mean, taken as a whole, becomes more and more symmetrical on either side 
of the centre of gravity. 

It may be remarked here, by the way, that a relation holds for the c. 
m. inequality, very similar to that which holds for the q. m. error, in a 
quantity X which is connected with other quantities x lf x 2 &c. thus, 

X = a 1 x 1 -\-a 2 x 2 -\-a & x s -\-&c, (61) 

where a 1} es 2 &c. may be essentially either + or — . The errors of x 1} x 2 
&c. are supposed to be independent, that is, the error of one has no influence 
on the error of another. If s 1} e 2 &c. denote the q. m. errors of x 1 , x 2 &c, 
and e denotes the q. m. error of X, then as is well known 

e 2 = (a l e 1 ) 2 +(« 2 s 2 ) 2 +Ke 3 ) 2 +&c. (62) 

(See the method of proof which I gave in Analyst, VIII, p. 139.) In 
like manner, it is easily seen that if the errors to which x x , x 2 &c. are liable 
are of such nature that + and — errors of equal amount are not equally 
probable, then denoting the c. m. inequalities of x u x 2 &c. by £ x , £ 2 &c> 
and that of X by £, we shall have by virtue of (37) 

C 3 = KCi) 3 +KC2) 3 +KC 3 ) a +&c (63) 

The a lf a 2 &g. are merely coefficients, so that any actual error of a x x x , 
for instance, is a x times the actual error of a^. Then the c. m. inequality 
of a x x t is OiCi, that of a 2 x 2 is a 2 £ 2 , and so on. The c. m. inequality of 
—x 2 being — £ 2 , that of (x t ±x 2 ) is the cube root of (£f ±Cl)- 

[To be continued.] 



